When did cosmic acceleration start? How fast was the transition? 
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Abstract 

Cosmic acceleration is investigated through a kink-like expression for the deceleration parameter (g). The new 
parametrization depends on the initial (qi) and final (g/) values of q, on the redshift of the transition from deceleration 
to acceleration (zt) and the width of such transition (r). We show that although supernovae (SN) observations 
{Goldl82 and SNLS data samples) indicate, at high confidence, that a transition occurred in the past {zt > 0) they do 
not, by themselves, impose strong constraints on the maximum value of zt. However, when we combine SN with the 
measurements of the ratio between the comoving distance to the last scattering surface and the SDSS+2dfGRS BAO 
distance scale (Sk/Dy) we obtain, at 95.4% confidence level, zt = 0.84±!3:i? and t = 0.51±!];?? for {Sk/Dy + Goldl82), 
and zt — 0.88±SJ;io and r = 0.35±o;lo for (Sk/Dv + SNLS), assuming qi = 0.5 and qj — —1. We also analyze the 
general case, g/ G (^co,0) finding the constraints that the combined tests {Sk/Dv + SNLS) impose on the present 
value of the deceleration parameter (go). 



1. Introduction 

Since the discovery of the accelerated expansion 
of the universe in 1998 |1I2) . considerable effort in 
cosmology has been devoted to determine the source 
of this acceleration. The two most common possi- 
bilities discussed in the literature are: the existence 
of an exotic component with sufficiently negative 
pressure (dark energy) and proper modifications of 
general relativity at cosmological scales (for recent 
reviews see [3]). 

One way of making progress in determining the 
cosmic expansion history is through a model by 
model analysis. Another is to carry out a phe- 
nomenological analysis with the use of different pa- 
rameterizations of the dark energy equation of state 
[4], the Hubble parameter [5] or the dark energy 
density [6j . This procedure may provide interesting 
pieces of information, but in general a parametriza- 
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tion assumes the existence of dark matter and dark 
energy as different substances (barring a few excep- 
tions no interaction in the dark sector is considered) 
and general relativity is in most cases assumed. In 
this framework, an important question regards the 
number of parameters necessary to get reliable con- 
clusions. If too many are used, the allowed region 
in the parameter space could be so large that it 
would not be possible to get firm conclusions [7]. 
Otherwise, if not enough parameters are used, the 
obtained results may be strongly dependent on the 
particular parametrization choice and misleading 
conclusions could be reached [8|. The strategy we 
follow here is to use a large (four) number of pa- 
rameters in order to be quite general, but, based on 
physical arguments, fix two of them from the start. 
We then relax the condition on one of the fixed pa- 
rameters and obtain the confidence surface on the 
other three. 

In this work we are mainly interested in the follow- 
ing questions: what is the redshift of the transition 
from decelerated to accelerated expansion? How fast 
was it? We investigate these by introducing a new 
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parametrization for the deceleration parameter (q) 
that depends on four parameters: the initial {qi) and 
final {qf) values of q, the redshift of the transition 
from deceleration to acceleration (zt) and a quan- 
tity related to the width in redshift of such transi- 
tion (t). With this formulation we aim to answer 
the above questions with the minimum amount of 
assumptions about the dark sector and the funda- 
mental gravitation theory. 

This paper is organized as follows: in Section 2, we 
present the new q parametrization and discuss some 
of its properties. In Section 3 the outcomes of the 
confrontation of this parametrization with two su- 
pernovae samples, the new Goldl82 SNe la from [9] 
and the first year data set of the Supernova Legacy 
Survey (SNLS) [lO], are obtained (first assuming 
g/ = — 1 and qi = 1/2). We show that current su- 
pernovae observations alone are not able to satis- 
factorily constrain the transition redshift. To break 
the SN la degeneracy, we combine this observable 
with the ratio of the comoving distance to the last 
scaterring surface {Sk{zis = 1098)) to the baryon 
acoustic oscillations (BAO) distance scale {Dy{z)) 
at zbao = 0.2 and zbao = 0.35 as estimated in 
[TT] . We show that the ACDM model is within the 
region allowed by SNLS+S'fc/I?^ results but is ex- 
cluded for Goldl82+Sk / data, at 95% confidence 
level. We then discuss the broader case with arbi- 
trary qf exhibiting the 95% confidence surface in the 
parameter space (zt,T,qo) {qg is the present value 
of q), obtained using the SNLS+ Sk/Dy data. Our 
conclusions are presented in Section 4. 

2. The Model 

At large scales, it is a good approximation to 
consider a spatially homogeneous and isotropic 
universe. With this assumption we are lead to the 
Friedman- Robertson- Walker metric: 



a{tf 
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1 - kr 



r^dn^ 



(1) 



where a{t) is the scale factor and k = — 1,0, +1 
characterizes the curvature of the spatial sections of 
space-time. From now on we will assume a fiat uni- 
verse (fc = 0), which is in agreement with CMB re- 
sults |T2] . 

In terms of the Hubble parameter {H = ^), the 
deceleration parameter can be written as: 
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Therefore, 

H — Hq exp 



(g(i)-M)dln(l + z) 



(3) 



In this work, we propose the following phe- 
nomenological functional dependence with redshift 
for the deceleration parameter: 

(g» " Qf) 
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where qi > (deceleration) and qf < (accelera- 
tion) are the initial (z ^ zt) and final {z = —1) val- 
ues of the deceleration parameter, respectively. The 
parameter zt denotes the redshift of the transition 
{q{zt) — 0) and r > is associated with the width of 
the transition. It is related to the derivative of g with 
respect to the redshift aX z = zt. More precisely. 
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T are demon- 



The infiuence of parameters 
strated in Fig. 

Expression (|4]) is similar in spirit to the one 
suggested in [13] (see also [7l8.|14j), but here we 
parametrize q{z) instead oiw{z). One of the advan- 
tages of using the above kink-like parametrization 
for the deceleration parameter is that zt has a very 
clear physical meaning. Different physical aspects 
and parameterizations of g(z) were also investigated 
in [15]. With the above definition, equation ^ is 
now integrated to give. 
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We now define an effective matter density parameter 
17,„_ EE lim ( '-^ I (1 + z)-2(i+90 _ (7) 



Ho 

where the limit should be understood as z >> zt- 

In most (and simplest) scenarios, in order to form 
large scale structures, the universe passes trough a 
kind of matter dominated phase such that, at early 
times (but after radiation domination), H^ oc (1 -|- 
z)^, which implies g = 1/2. In this work we fix qi = 
1/2 reducing to three the number of free parame- 
ters. In principle, with this assumption we are losing 
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Fig. 1. Influence of parameters zt and r in the functional form of the deceleration parameter for the special case qi = 0.5 
and qj = —1. Left - q{z) for zt = 1.0 and t=0.1 (green dashed), 0.3 (full black) and 0.5 (blue dot-dashed). Right - q{z) for 
T = 0.3 and zt = 0.5 (green dashed), 1.0 (full black) and 1.5 (blue dot-dashed). 



generality but the question is: how much qi can devi- 
ate from 1/2 during large scale structure formation? 
In the general relativity framework, in some models 
with a constant coupling {5) between dark matter 
and dark energy this condition {qi = 1/2) is not sat- 
isfied. In this case we have H'^ oc (l-f during 
matter domination. But, what are the allowed val- 
ues for 51 In [16] it has been shown that background 
cosmological tests, impose \5\ < 0.1. Taking into ac- 
count matter perturbations stronger constraints on 
the coupling can be obtained [iTj. Another possibil- 
ity is to consider models in which matter has pres- 
sure, such that when it dominates qi ^ 1/2. How- 
ever, if matter perturbations are adiabatic, due to a 
finite speed of sound, the mass power spectrum will 
present instabilities ruling out these models unless 
p = {q = 1/2) or very close to it. In principle, it is 
possible to circumvent this kind of problem by as- 
suming entropy perturbations such that Sp = [18] . 
However, in this case the models may have prob- 
lems with lensing skewness as pointed out in [19] . 
Although there are some indications that by fixing 
qi = 1/2 we are not losing much in our description of 
the majority of the viable models, relaxing this con- 
dition requires further investigation, and we leave it 
for future work. 

In the specific case qi — 1/2 we have , 

/ -I .-r{l-2qf) 

ri^oo- (^1-^(1 + Zt)'^"J • (8) 

With the above definition, we can eliminate Zt from 
equation ^ and rewrite it as 
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X f anoo"^'^' + (1 - n^'^''''){l + z)-i] (9) 

The above expression for H{z) in terms of ^moo is 
very useful to make connections with models already 
discussed in the literature. For instance, it is simple 
to verify from ^ that the parametrization (j4]), in 
the special case qi ~ 1/2, is related to the "Modi- 
fied Polytropic Cardassian" (MFC) model [20]. This 
model depends on three parameters: m (denoted by 
q in [20]), n and ^mo- If we identify, ^mo = f^moo, 
m = 1/(t(1 - 2qf)) and n = 2/3(1 + g/), it fol- 
lows that the two models have the same kinemat- 
ics. Note that, since g/ < 0, the condition n < 
2/3 follows naturally. We remark that in the MFC 
model, Qrno is the present value of the matter den- 
sity parameter, while flmoo is defined at high red- 
shift. These two quantities do not necessarily have 
the same value in the general case if, for instance, 
dark matter and dark energy are coupled [21]. As an 
example, consider models with a variable coupling 
between dark matter and dark energy (assumed to 
have constant equation of state Wx), and such that 
Px/Pm = Pxq/ PmoO'^ [22]. These models can be de- 
scribed by Eqn. ^ if we identify flmo = ^Inoc^ ^'^^ , 
^ = l/r and wx = —(1 — 2g/)/3. As remarked be- 
fore, in our formulation it is not necessary to make 
strong assumptions about the dark sector or gravity 
theory. In MFC model the universe components are 
specified to be matter and radiation; there is no dark 
energy. The parametrization (|4]) includes the MFC 
model (and the coupling models above) as special 



cases. 

Neglecting baryons, the quartessence Chaplygin 
model {p = -M-*("+i)/p") is obtained if we as- 
sume Qi = 1/2, Qf — —1, identify 1/t ~ 3{l + a) and 
^mo. = (1 - i«o)'/^'+"\ where Wo = -MVpo+' is 
the present value of the equation of state parameter. 

The conventional dark energy model with con- 
stant equation of state {wx) is obtained if we iden- 
tify flmoo = f^mo and impose the condition —3wx = 
1/t = (1 — 2(7/) in Eqn. In particular, if g/ — —1 
and T — 1/3, ACDM is recovered. For this model the 
transition redshiftis equal to {2{l—flmo) / ^moV^^ — 
1. Identifying ACDM in the parameter space is very 
convenient; it fits current data quite well and we 
should expect the "true" cosmology not to be far 
from this limit. We remark that, in the framework of 
general relativity with non-interacting dark matter 
and dark energy, if r < 1/3 and g/ — —1, the dark 
energy component will present a transient phantom 
{w < —1) behavior, that could either have started 
in the past or in the future (z < 0). Models with 
T > 1/3 are always non-phantom. 

It is curious that if we apply the definition of t, 
given by Eqn. ([5]) (assuming qi = 1/2 and g/ = —1), 
to the flat DGP brane-world model [24J we obtain 
r = 1/2, independent of ilmo- Therefore q-models 
with zt = (2(1 - nmo)V ^m^Y'^ - 1 (the DGP 
redshift transition) and r « 1/2 are expected to be 
a good approximation for flat DGP models. 

3. Observational Constraints 

One of the main questions today in cosmology is 
to know if cosmic acceleration is generated by a cos- 
mological constant or not. The data seem to indicate 
that models "close" to ACDM are favored. In our 
analysis we flrst consider the special case of models 
that have a flnal de Sitter phase (g/ — —1). In this 
case, the flat ACDM model is more easily identifled 
in the parameter space allowing a simple test of the 
ACDM paradigm. The more general case, with ar- 
bitrary g/, will also be briefly considered. 

Assuming g/ = — 1 we now derive constraints on 
the parameters r and zt by combining supernovae 
measurements with the ratio of the comoving dis- 
tance to the last scaterring surface, Sk{zis = 1098), 
to the BAO distance scale, Dy{z), at zbao = 0.2 
and Zbao — 0.35, as estimated in [11]. In fact, the 
ratio Sk/Dy times zbao is equal to the ratio of the 
CMB shift parameter (7^) [26J at z^ to the BAO pa- 
rameter A{zbao) [25] ■ This observable is appropri- 



ate for our purpose for two reasons. First, it does not 
explicitly depend on the exotic dark constituents of 
the universe and neither on the gravity theory. It 
is essentially controlled by the function H{z)/Hq. 
Second, complementarity with the SN constraints 
is generated because the Sk/Dy ratio and SN are 
sensitive to distances to objects (events) in different 
redshift range; with supernovae we are measuring 
distances up to z ~ 1 — 2, while Sk depends on the 
comoving distance to z 1100. 

In the SN la analysis we considered both, the 
Goldl82 [9] and the SNLS tlOj samples. To deter- 
mine the likelihood of the parameters we follow the 
same procedure described in these two references. 
In our computations, when marginalizing over the 
Hubble parameter, we use a Gaussian prior such that 
h = 0.72 ± 0.08 [27]. In Figl2] (left-panel), assuming 
g/ — —1, we display constant confldence contours 
(68% and 95%) in the (arctanzt, t) plan allowed by 
SN experiments. Notice that, for the two SN data 
sets Zt < is not allowed at a high confldence level, 
indicating that a transition occurred in the past. We 
remark that this is expected in ACDM models (or 
other models) that have a non-null transition time, 
but our results indicate that this is true even if the 
transition is instantaneous (r = 0) . This conclusion 
also appHes if g/ ^ —1. Furthermore, it is also clear 
in Figl2] (left-panel) that current SN observations 
cannot impose strong constraints on the maximum 
value allowed for zt . Since SN observations prove the 
universe only up to redshift z ~ 1 — 2, in a model 
in which the transition is slow (r > 1), even if zt is 
high, the distance to an object, let say, at z < 1, can 
be similar to the distance to the same object in an- 
other model in which zj < 1 with a faster transition 
(smaller r). This explains the shape of the SN con- 
tours. By comparing the confldence contours for the 
two data sets, we observe that those from Goldl82 
are shifted to lower zt with respect to those from 
SNLS. We remark that even in the region of more 
interest (zj < 1), the difference between the out- 
come of the two SN la data sets, although not so 
severe, exists and is important. Similar results were 
obtained in |28|, which can be related to possible in- 
homogeneities present in the Goldl82 sample and 
should be further investigated. 

To obtain the constraints on the parameters from 
the Sk/Dy test, we use a statistics taking into 
account the correlation matrix and the ratio Vg/Sk 
given in [TTj. Since we are assuming flat space we 
have, Dy{zBAo) = [zbaoH-\!;^''" dzHfY'^ 
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Fig. 2. Left - Constraints imposed by SN la observations. The contours represent 68% and 95% c.l. . The green-dashed curves 
are from Goldl82 data set and the blue-solid ones stand for SNLS first year data set. Right - Constraints imposed by S^/Dv 
measurements (68% and 95% c.l.) as explained in the text. The red horizontal line in both panels corresponds to flat 
ACDM models. Notice that, for these models, higher values of zt correspond to smaller values of OmO. For both panels qf = —1 
is assumed. 



and Sk = /q^°^^ dzH~^. The He term in the defini- 
tion of Sk incorporates the necessity of taking into 
account the contribution of radiation at very early 
times, which is not included in our parametriza- 
tion Eqn. (0]). One may argue that by introducing 
radiation we are losing generality. However, any 
late time modification of the standard cosmologi- 
cal model should satisfy Big-Bang nucleosynthesis 
(BBN) constraints at very early time. We know 
that radiation exists and what should be the de- 
pendence of the Hubble parameter with redshift 
at very early times (when radiation dominates) in 
order not to spoil BBN's success. During this phase 
oc (1 -t- z)"^ and q = I. Therefore, to take radia- 
tion into account, we add the term flroi^ + to 
the right hand side of ^ when applying it to cal- 
culate Sk - If we do not consider it, we would have a 
~ 18% error in estimating Sk- 

We marginalize the likelihood over h with the 
same Gaussian prior used in the supernovae analy- 
sis. In fact, Sk is almost independent of h; the de- 
pendence entering only through the radiation term. 
In FigO (right-panel) we show constant confidence 
contours (68% and 95%) in the {zt,T) plan allowed 
by the Sk/Dy test. It is worth to be mentioned that 
the shape of the 95% contour is similar to contours 
of constant flmoo and we can think that this test 
essentially constrains this quantity. The same kind 



of behavior also appears in fiat, constant w mod- 
els [H]. Furthermore, it is clear from the figure the 
complementarity between the SN and this test. 

To get the combined (SN+Sk/Dy) results we mul- 
tiply the marginalized likelihood functions. In FiglS] 
(left-panel) we show the results (68% and 95% c.l.) 
of the Sk/Dy test with the Goldl82 (green dashed 
contours) and with the SNLS (blue dot-dashed con- 
tours) data set. The red horizontal line in the figure 
represents the ACDM limit (r = 1/3). It shows that 
this model is in good agreement with the Sk / Dy + 
SNLS data. After marginalizing over the extra pa- 
rameter L29J we find for Sk/Dy + Goldl82 (at 95.4% 
confidence level), Zt — 0.84±g;}7 and r 
while for S'^/D^+SNLS we have, zt 
T — 0.35±q;Jo- For a model with q/ = —1, r = 
0.35 and zt — 0.88, we obtain from Eqn. ([8]) that 
^moo — 0.23. It is also simple to show that the age 
of the universe in this particular model (assuming 
h — 0.72) is 14.0 Gyr and that cosmic acceleration 
started 7.2 Gyr ago. Notice that the ACDM cos- 
mology (r = 1/3) is in good agreement with Sk/Dy 
+ SNLS, but excluded at 95.4% confidence level by 
the Sk/Dy+Goldl82 data. This discrepancy reveals 
tension between the two SN data samples and re- 
inforces the necessity of better SN data to clarify 
the issue. We also display in the same figure (solid 
contour) what should be expected from future sur- 



0.51±g;?f, 
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Fig. 3. Left - 68% and 95% confidence levels imposed by the combined data sets of figure JSJ . The green-dashed (blue-dot-dashed) 
contours represent Goldl82 (SNLS)+ S^/Dy. The small black-solid contour (95% c.l.) was obtained from simulated data as 
explained in text. For the figure Qf = —1 is assumed. Right - We show, in the parameter space {zt,T,qo), the 95% confidence 
surface for the qf general case (qj g (— oo,0)), obtained using SNLS+5fc/Di, data. 



veys when combining SN la+Sk/Dy measurements. 
In our Monte Carlo simulations we used as fiducial 
model a flat ACDM model with = 0.23 (r = 
1/3, zt ~ 0.88). For SN la we considered a SNAP- 
like survey assuming that the intercept is known. For 
the Sk/Dy test we used a conservative (but some- 
what arbitrary) assumption that the uncertainties 
will be reduced to 2/3 of their current values. We 
also assumed that the correlation coefficient would 
remain the same. In the figure we show the 95% con- 
fidence contour. 

We also analyzed the broader case with arbitrary 
qf. Our parametrization ^ allows us to determine 
the present value of q (go) in terms of zt, r and qf. 
Although the considered data sets do not impose a 
lower bound for qf, they do constrain go (we found 
— 1.4 < go ^ —0.3). In Figl3] (right-panel) we show, 
in the parameter space {zt , t, go) , the 95% confidence 
surface for the general case (g/ G (—oo, 0)), obtained 
using SNLS+Sk/Dy data. 



models like fiat ACDM, DGP and others. We con- 
firmed that there is a tension between Goldl82 and 
SNLS data sets with a quite general formulation. We 
also exhibited what should be expected from future 
SN + Sk/Dy observations and, relaxing the condi- 
tion qf — —1, obtained current constraints on the 
parameters zt, r and go. A more detailed analysis 
of the consequences of our parametrization in the 
case g/ ^ —1 is still necessary. For instance, our 
parametrization is not able to describe, in all their 
redshift range, models that are now accelerating but 
that decelerates again in the future [30] . The descrip- 
tion of the expansion history of these models is more 
complicated since it requires a second transition. It 
would be interesting to investigate under what con- 
ditions it would be possible to describe, with our 
parametrization, the behavior of these kind of mod- 
els for z > 0. The case qi ^ 1/2 should also be fur- 
ther investigated. These issues will be discussed in 
subsequent work. 



4. Conclusion 



In this work, with a formulation that avoids strong 
assumptions about the dark sector and / or the met- 
ric theory of gravity, we showed that by using only 
SN data the transition redshift (from decelerated 
to accelerated expansion) could be very large {zt > 
10). We demonstrated the importance of combin- 
ing the SN test with the Sk/Dy test to better con- 
strain the parameters zt and t. We introduced the 
parameter r and showed its relevance to characterize 
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